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Abstract 



Qj \ A study of the chiral behavior of pseudoscalar masses and decay constants is 



carried out in quenched lattice QCD with Wilson fermions. Using the modified 



^ ' quenched approximation (MQA) to cure the exceptional configuration problem, 



accurate results are obtained for pion masses as low as ~ 200 MeV. The anoma- 
lous chiral log effect associated with quenched r/' loops is studied in both the 
relation between vs. niq and in the light-mass behavior of the pseudoscalar 
and axial vector matrix elements. The size of these effects agrees quantitatively 
with a direct measurement of the rj' hairpin graph, as well as with a measure- 
ment of the topological susceptibility, thus providing several independent and 
quantitatively consistent determinations of the quenched chiral log parameter 
6. For j3 = 5.7 with clover-improved fermions {Csw = 1-57) all results are 
consistent with 6 = 0.065 ± 0.013 . 



1 Introduction 



Until recently, efforts to study the chiral limit of lattice QCD with Wilson-Dirac fermions 
have been impeded by the large non-gaussian fluctuations encountered at light quark mass 
(the "exceptional configuration" problem). These statistical problems have been shown to 
be symptomatic of the nonconvergence of the quenched lattice Monte Carlo integration due 
to the presence of exactly real Dirac eigenmodes in the physical mass region [Q]. These 
real eigenmodes are of topological origin but are displaced from zero mass by the explicit 
chiral symmetry breaking of the Wilson-Dirac operator. The modified quenched approxi- 
mation (MQA) is a prescription for constructing a modified quark propagator by shifting 
positive mass real poles to zero mass in a compensated way.[Q] This procedure is the most 
straightforward and effective way of removing the lattice artifact which produces excep- 
tional configurations while leaving the results for heavier quarks outside the pole region 
essentially unchanged. In the continuum limit, the spread of eigenmodes to positive mass 
shrinks to zero, so that if we shift all poles above a given quark mass, the fraction of gauge 
configurations in an ensemble that would require MQA pole-shifting approaches zero in 
the continuum limit. In this sense, the MQA has the right continuum limit, and provides 
a sensible definition of the otherwise undefined quenched theory at finite lattice spacing 
and small quark mass. Exactly how effective the MQA is in restoring the chiral behavior 
expected from quenched continuum arguments is a question that can best be addressed by 
detailed calculations. As we will show in this paper, by using the MQA procedure we are 
able to study the behavior of pseudoscalar masses and decay constants with high statistics 
down to much lower pion mass (ca. 200 MeV) than has previously been accessible at /? = 5.7 
(for example, a recent high statistics study of quenched chiral behavior extends only to 
~ 300 MeV). The results confirm in considerable detail the expected chiral behavior 
predicted for these quantities from chiral Lagrangian arguments. The overall size of the 
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quenched chiral log parameter 6 is about a factor of two smaller than the 5 ~ 0.17 expected 
from the theoretical result 

5= (1) 

using the physical values of m^i and f-,^. (Here the normalization of /^r is such that its 
physical value is ~ 95 MeV. Chiral perturbation theory and the physical mass of the rj' 
give rriQ ~ 850 MeV.) Nevertheless, the quantitative relations between various quantities 
inferred from continuum chiral Lagrangian arguments appear to be well-satisfied. In fact, 
the smallness of 5 is consistent with our previous calculations at heavier quark mass which 
did not use MQA shifting [|]. 

In this paper, we describe several independent quantitative estimates of the quenched 
chiral log (QCL) parameter 6: (1) a direct measurement of the r/' hairpin diagram, (2) a 
calculation of the topological susceptibility combined with the Witten-Veneziano formula 

(3) measurement of the QCL effect in vs. ruq and in the vacuum-to-one-particle 
matrix element of the pseudoscalar density, (4) by fitting cross-ratios of masses and matrix 
elements using the lowest order chiral perturbation theory prediction (to first order in 5), 
and (5) by a global fit of all masses and matrix elements to the next-to-leading order 
chiral perturbation theory prediction. When the two quark masses in the pion are allowed 
to vary independently, as in (4) and (5), the corresponding masses and decay constants 
provide additional tests of quenched chiral log structure (c.f. the analysis of the CPPACS 
collaboration||3|). In this paper we present the results of analyses of chiral log structure both 
in the "diagonal" (equal quark mass) sector, and for mesons with unequal quark masses. 
The various determinations of 5 are all consistent within statistics, giving an exponent 
6 = .065(13) for clover-improved Wilson fermion action, and 6 = .06(2) for unimproved 
Wilson. 

As discussed below, operators which do not change chirality, such as the axial-vector 
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current, are not expected to exhibit any QCL singularity (for equal mass quarks). In 
agreement with this theoretical expectation, our results for the mass dependence of the 
axial-vector decay constant show smooth analytic behavior with no enhancement in the 
chiral limit. The contrast between the singular chiral behavior of the pseudoscalar matrix 
element and the smooth analytic behavior of the axial- vector matrix element is a particularly 
convincing piece of evidence that we are indeed seeing the effects of virtual rl loops. 

In Section 2 we briefly review the essential features of quenched chiral behavior in QCD 
in the language of effective field theory. Details of the lattice calculations are presented 
in Section 3. The direct determinations of b from the ry' hairpin correlator and from the 
topological susceptibility are presented in Section 4. In Section 5, QCL effects in the 
squared pseudoscalar mass as a function of quark mass are discussed, for "diagonal" mesons 
composed of equal mass quarks only. The QCL behavior of the pseudoscalar and axial- 
vector decay constants(again, for diagonal mesons only) is studied in Section 6. In Section 

7 we extract the chiral log parameter b from an analysis of cross-ratios of pseudoscalar 
meson masses and decay constants, along the lines of the CPPACS analysis |p. In Section 

8 we compare our lattice data for all mesons (equal and unequal quark masses) with the 
results of a detailed quenched chiral perturbation theory (to order -p^) calculation for the 
pseudoscalar masses and decay constants. Finally, in Section 9 we summarize our results 
and draw some broad conclusions. Some technical remarks on the efficacy and accuracy of 
the all-source technique used to extract hairpin amplitudes are relegated to Appendix 1, 
while chiral perturbation theory formulas for pseudoscalar and axial decay constants are 
given in Appendix 2. 
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2 Brief review of quenched chiral theory 

The presence of anomalous chiral behavior induced by the quenched approximation was 
first investigated in the work of Sharpe and Bernard and Goltermann Q . Physically, the 
quenched chiral log effect can be understood as a consequence of the absence of topological 
screening in quenched QCD. Diagramatically, the quenched approximation discards some, 
but not all, of the hairpin mass insertions. Instead of the massive r/' propagator of full 
QCD, the quenched hairpin diagram exhibits a double Goldstone pole, see Fig[j^. This 
leads to infrared singular r( loop diagrams, such as that shown in FigQ. which alter the 
chiral behavior of certain matrix elements. A convenient way to understand the effect of 
the QCL singularities is to interpret them as a renormalization factor in the chiral field. 
Begin with a C/(3) x U (3) chiral Lagrangian with a chiral field 



V = exp 



L i=0 



where Aq = 1 and ttq = rj' represents the SU(3)-flavor singlet meson. Now consider the 
effect of integrating out the r]' field. The remaining SU{'i) x SU (3) Goldstone fields will be 
renormalized 

U ^ (e-o//)c7 = exp [-(^o')/2/'] U (3) 
where, in the quenched theory, 



d^p 1 f d^p 



^ ^ J {2^Y p2 + ^2 ^ ] (2vr^4 (^2 + ^2)2 (4) 

In full QCD, the r}' mass is large, so the effect of integrating out the r]' field is to induce 
a finite renormalization that is nonsingular in the chiral limit. The chiral limit is described 
by the S'[/(3) x SU{3) chiral Lagrangian without an r]' . Now consider the effect of quenching 
on the result of integrating out the r]' field. The double Goldstone pole in the quenched ij' 
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= X 

Figure 1: Chiral Lagrangian representation of the quark hairpin diagram as an 77' 
mass insertion. 




Figure 2: A quenched 77' loop diagram which produces anomalous chiral log behavior. 
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propagator produces logarithmically divergent loop integrals in cliiral perturbation theory. 
This results in a renormalization factor for the field U that is singular in the chiral limit, 
and thus alters the m,r behavior of the chiral field: 

U exp [-,51ogm^l U = U (5) 

where U is nonsingular as m-Tr 0. In the following subsections, we will discuss the chiral 
log effect (or lack thereof) in the pseudoscalar and axial- vector charge matrix elements and 
in the mass of the pion as a function of quark mass. 
Define the decay constants fp and Ja 

(0|Vi7'V'k(p)) = fp 

{0\i^^'>ri^\7rip)) = p^fA (6) 

Prom the chiral-field expressions for the quark bilinears, 

ipj^tp oc U-U^ 
V^tVV' oc i(u-'^di'U - {d''U-^)u) , (7) 

we conclude that a singular chiral log factor will appear in the quenched calculation of fp, 
but not of fA, 

/•quenched f ^ \ f 

/•quenched I /on 

Ja = M [o) 

where fp and fA go to a constant in the chiral limit. Both of these expectations are 
confirmed by the data, as discussed below. The chiral behavior of the pion mass as a 
function of quark mass is easily derived from the results for the pseudoscalar and axial- 
vector densities, combined with PCAC. This gives 

mi ~ const. X mn X ( — ^r)^, m„ ^ (9) 
mi 



This predicted behavior is also confirmed by the lattice results, as discussed in Section 5. 

The quantitative significance of the quenched chiral log behavior that we observe is 
further reinforced by comparing it with a direct calculation of the anomalous exponent S. 
We have done this calculation in two related but distinct ways. One is a calculation of the 
hairpin diagram Fig[|^. 

Ah{x) = {Tt^^G{x, x)Tn^G{0, 0)) (10) 

where G{x, y) is the quark propagator. The size of the hairpin mass insertion determines 
the coefficient of the logarithmic divergence in the ry' loop diagrams. Let us denote by 
mg the value of the r/' mass insertion vertex obtained from the long-range behavior of the 
hairpin diagram. Then the predicted chiral log exponent is given by Equation (|l|). 

The other quantity which provides an evaluation of 6 is the topological susceptibility, 
Xt- The Witten-Veneziano formula allows us to relate xt to the r/' mass insertion, 

ml = -jfxt (11) 

where Nf is the number of (light) quark fiavors. The topological susceptibility is given in 
terms of winding numbers on a lattice of four- volume V by 

..^^ (12) 

To calculate winding numbers, we use the integrated anomaly technique ("fermionic 
method") of Smit and Vink The configuration- by-configuration calculation of wind- 
ing numbers provides a determination of the xti and hence via (11) and (|^), an estimate of 



6. The results obtained from the hairpin residue and from the topological susceptibility are 
in good agreement for clover improved fermions, as discussed in Section 4. 
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3 Lattice calculation of pseudoscalar masses and 
decay constants 

Most of the results which we will focus on in the subsequent discussion are obtained from the 
"b" ensemble of gauge configurations from the Fermilab ACPMAPS library. This ensemble 
consists of a set of 300 quenched configurations on a 12^ x 24 lattice at /3 = 5.7. The 
fermion action used in our calculations was clover improved with a clover coefficient of 
Csw = 1-57. To get some estimate of lattice spacing effects and finite volume corrections, we 
also calculated masses and decay constants on 200 configurations of the b-lattice ensemble 
with unimproved Wilson fermions {Cgw = 0), and on the a-lattice ensemble (16^ x 32 at 
(3 = 5.7), also with unimproved Wilson fermions. The masses and decay constants and their 
quoted errors are obtained from fully correlated fits to the smeared-smeared and smeared- 
local propagators, using a smeared pseudoscalar source in Coulomb gauge combined with 
local pseudoscalar and axial-vector sources. Calculations with different smearing functions 
give results that are consistent within less than one standard deviation, indicating that 
the systematic error associated with excited state contamination is less than our statistical 
errors. For the hairpin correlator, the issue of excited state contamination is addressed 
in detail in Section 3.2 by comparing the ratio of smeared and local hairpin propagators. 
Remarkably, we find that for small quark masses there is very little excited state contribution 
to the hairpin propagator, even at time separations as small as t = 1. This is in marked 
contrast to the valence pion propagator, which has substantial excited state components at 
short time. We conclude that the hairpin vertex itself is largely decoupled from the excited 
states in the pseudoscalar channel. 

The essential new ingredient in the present study which allows a detailed exploration 
of chiral behavior is the use of the pole-shifting procedure of the modified quenched ap- 
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proximation. The details of this procedure and its effectiveness in resolving the exceptional 
configuration problem have been described previously For each gauge ensemble and 
choice of fermion action, we carried out a careful scan of each configuration for poles over 
a range of quark masses, starting at a value heavy enough to be beyond all real eigen- 
mode poles and going to nearly zero mass. For example, in the run with clover action, 
where Kc = 0.14329, poles appeared as low as k ~ 0.1417 (m^ ~ 45 MeV), and we scanned 
for and located all poles up to 0.1431 . The value of the integrated pseudoscalar charge 
Qb = J d'^xip'y^ip is calculated for a sequence of hopping parameter values, using the same 
allsource method that is used to calculate the hairpin propagator (see below). A Fade fit to 
these values determines the location of any poles within and somewhat beyond the range 
scanned. Extremely precise pole locations can be obtained by performing further conjugate 
gradient inversions very close to the pole positions determined by the Fade fit. Using the 
stabilized biconjugate gradient algorithm, we are able to perform inversions at hopping pa- 
rameters very close to the pole position without any major increase in convergence time. In 
our calculations, we have located the pole positions as a function of hopping parameter k, to 
at least 8-digit accuracy in all cases. Once all the visible poles in an ensemble are located, 
their residues in the quark propagator are determined by performing inversions slightly 
above and below the pole and subtracting. We found that accurate pole residues could be 
computed for a pole at kq by inverting at k = ± .000001. (The computation of pole 
residues can be done very economically by noting that the pole contribution to all 12 color- 
spin components of the quark propagator can be obtained from a single color-spin inversion 
above and below the pole, i.e. only 1/12 of a full propagator calculation is required.) 

An alternative procedure for locating poles of the quark propagator oc l/(i/* — m^) is to 
use the Arnoldi algorithm Q to partially diagonalize the Wilson-Dirac operator in the region 
around zero mass. The returned eigenvalues are the pole positions, and the residues needed 
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to perform the pole-shifting procedure may be reconstructed from the Arnoldi eigenvectors. 
This method locates not only the real eigenvalues but also the complex ones in the continuum 
band. The Arnoldi analysis has been carried out on a subset of the configurations used in 
this investigation, extracting approximately 50 eigenvalues from each gauge configuration 
analyzed. For the lattices which exhibited visible poles in the bicongradstab scanning 
procedure described above, the results of the Arnoldi calculation agree accurately with 
that analysis for both the pole positions and residues. The low-lying spectrum for an 
"exceptional" gauge configuration {b : 021000) is shown in FigQ. (Note: The positive mass 
region is to the right in this plot. The three modes farthest to the right on the real axis 
were MQA shifted.) 

The effect of the MQA pole-shifting procedure is to eliminate the problem of exceptional 
configurations and to dramatically improve error bars on all quantities calculated from 
light quark propagators, as shown in detail in Ref.|||]. Since the rj' hairpin propagator is 
particularly sensitive to the topological structure of gauge configurations, the improvement 
obtained by using the MQA for the hairpin calculations is even more striking than that for 
valence-quark meson propagators. The MQA-improved results are sufficiently accurate to 
allow a detailed study of the time-dependence of the hairpin propagator even as far out as 
t = 9 or 10. Over the entire time range starting from t = 1, this time-dependence is in 
reasonably good agreement with that predicted from a pure double pion pole form. The lack 
of any significant single-pole contribution to the hairpin propagator is consistent with the 
absence of excited states. (A heavy excited state on one side of the hairpin vertex would 
produce an effective single pion pole term.) A single-pole term in the time-dependence 
could also arise from a p'^-dependent hairpin vertex. In the chiral Lagrangian, this would 
correspond to a renormalization of the r/' kinetic term in addition to a mass insertion. The 
fact that the time-dependence of the hairpin correlator is well-described by the pure double 
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Figure 3: A portion of the low-lying spectrum of the Cgw = Wilson-Dirac spectrum 
of a quenched gauge configuration (b_021000 of the ACPMAPS library) in the region 
A ~ 1/ Kc- This configuration has four real poles in the continuum band. The vertical 
line is at i?e A = 1/nc- The three real poles farthest to the right on the plot require 
shifting. 
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pion pole form over the entire time range from i = 1 to 12 lends quantitative support to the 
assumption, often used in phenomenological discussions, that the momentum dependence 
is mild and the r/' hairpin can therefore be treated simply as a mass insertion. The time- 
dependence is discussed in more detail in Section 4.2. The final results quoted for the mass 
insertion value rriQ are obtained from a 1-parameter fit of the hairpin to a pure double-pole 
formula, with the pion mass held fixed at the value obtained from the corresponding valence 
pion propagator. 
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4 The rj mass and the chiral log parameter b 



4.1 Topological susceptibility 

For both the hairpin correlator and the calculation of the integrated pseudoscalar charge 



Qs, the method used is one introduced into such loop calculations in Ref.|ll|. This method 
employs an "allsource" quark propagator calculated with a source that consists of a color- 
spin unit vector on all sites of the lattice. This allows closed quark loops originating 
from any space-time point to be included in a calculation (e.g. of = TrG^^ or of a 
hairpin diagram), relying on random gauge phases to cancel out the gauge- variant open 
loops. Even on a single gauge configuration, this method is a reasonably accurate way of 
calculating global quantities like Q5, since random phase cancellation of noninvariant terms 
should take place in the sum over sites. The topological winding number u of each gauge 
configuration can be determined using the integrated anomaly equation [^, 

1/ = -imqQ^ (13) 

Thus we expect the quark-mass dependence of Q5 to exhibit a simple pole at = with 
residue given by the winding number. The plots of Q5 vs. niq for two typical configurations 
(after MQA pole shifting) are shown in Fig[p. The solid line in each case is the best 
single-pole fit. For the 300 lattice ensemble with clover action, the distribution of winding 
numbers determined in this way is shown in Fig[^. 

The topological susceptibility xt can now be calculated as the mean squared winding 
number per unit four-volume 

X, = ^ = ^(Q?) (14) 
By calculating the last expression in (14) we observe only a slight quark- mass dependence of 
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Figure 4: as a function of quark mass for two typical gauge configurations after 
MQA pole shifting. 
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Figure 5: The distribution of winding numbers determined from the integrated pseu- 
doscalar density. 

the result, as shown in Fig|Q. Extrapolating to the chiral limit, we obtain, for Csw = 1-57, 

Xt = 6.48(58) X IQ-"^ = (188 MeV)^ (15) 

where the first result is in lattice units, and the second is obtained by using the charmonium 
scale at /3 = 5.7 of = 1.18 GeV. The corresponding result for unimproved Wilson 
fermions {Csw = 0) is 

Xt = 3.24(39) X 10"^ = (158 MeV)^ (16) 

Using the topological susceptibility and the value of axial vector decay constant obtained 
from the valence propagator fits (see Section 6), the Witten-Veneziano formula gives 

5 = .063(6) (17) 
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Figure 6: Topological susceptibility calculated by the fermionic method on a 12^ x 24 
lattice with clover improved {Csw — 1-57) fermions as a function of the quark mass 
used to determine Q^. 
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for clover improved quarks {Csw = 1-57), and 

6 = .074(9) (18) 



for unimproved Wilson quarks. In ( [17|) and (|1^) we used the bare quark mass obtained from 
the hopping parameter to determine winding numbers from values. Here and elsewhere, 
we have taken the bare quark mass to be the pole mass 

m, = log (^1 + ^(^-1-^-1)) (19) 

but the results for the chiral log parameter are not significantly different if we use the naive 
bare mass {k^^ — k'^^)/2. If instead we use the current algebra mass 

m^^ = fAml/2fp , (20) 

the result (|l7|) is essentially unchanged, while the Csw = result (18) is decreased to 

6 = .032(4) (21) 

(See discussion following Eq. (|28|).) 

An important check on the calculation of topological susceptibility is to show that (i/^) 
has the proper dependence on volume, i.e. {u"^) oc V. In Fig0 we compare the topological 
susceptibility calculated on a 12'^ x 24 lattice and on a 16^ x 32 lattice, both with unimproved 
Wilson fermions {Cgw = 0). The results shown in FigQ are consistent with the expected 
extensive property, i.e. a linear volume dependence. 

4.2 Hairpin correlator and the r]' mass insertion 

Using the same allsource propagators as in the previous subsection, we calculate the hair- 
pin contribution to the flavor singlet pseudoscalar propagator, i.e. the loop-loop correlator 
(p!o|). Earlier calculations of this correlator ^, 12 1 were restricted to relatively heavy 
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Figure 7: The volume dependence of (z/^) for 12"^ x 24 and 16^ x 32 lattices, for 
Wilson-Dirac fermions, Cg^ = 

quark mass and had large errors which prevented a detailed study of time-dependence. As 
discussed in the Introduction, the statistical problems encountered in these previous inves- 
tigations arise from exactly real Wilson-Dirac eigenmodes, the effect of which is magnified 
by the fact that the hairpin propagator receives its largest contributions from topologically 
nontrivial gauge configurations which necessarily contain such real modes. The MQA pole- 
shifting procedure is thus particularly effective in improving the hairpin calculation. Fig|p 
shows an example of a hairpin propagator before and after MQA improvement of the cor- 
responding quark propagators. The quark mass is still rather heavy here (m^ ~ 36 MeV, 
rnvr = .386a^^ ~ 450 MeV). For even lighter quarks, the unimproved hairpin propagator 
is unmeasurable, with errors > 100%, while the MQA improved hairpin is still quite well 
determined, allowing a reasonably accurate measurement of the ij' mass insertion even at 
the lightest quark mass we have studied. 
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The size and time-dependence of the hairpin correlator is measured accurately enough in 
the MQA method to address the two issues of excited state contamination and p^-dependent 
vertex insertion terms mentioned in the Introduction. These two effects are distinct, but 
they are difficult to disentangle from the time-dependence alone, since they both have the 
effect of adding a single pole term to the correlator. Fortunately, there is another way to 
determine the presence or absence of excited states, namely, to study the ratio of hairpin 
correlators obtained from smeared and local V'7^V sources. This can be compared with the 
overlap of the same smeared and local sources with the ground-state pion, as determined 
from the laige-t behavior of the corresponding valence pion propagators. 

We construct smeared-source hairpin correlators by a modification of the allsource 
method used for the local-source hairpins. In the latter, the source used for propagator 
inversion was a unit color-spin vector on every site. In order to obtain meaningful results 
for smeared source propagators, we must perform the smearing in Coulomb gauge. The 
smeared sources used for the valence pion propagators were constructed using an exponen- 
tial smearing function oc e^'^^' . Based on other studies of hadron wave functions at this 
value of P, we took A = 0.5 in lattice units. There is an additional subtlety in the imple- 
mentation of Coulomb gauge smearing in the allsource method. Since this method relies on 
random gauge phase cancellations, the actual sums over sites for the two ends of the hairpin 
must be carried out in the original unfixed gauge. In fact we carry out all calculations in 
the unfixed gauge, just as in the local calculation. The only difference is that the source 
used for propagator inversion is a "smeared allsource" which is constructed by the following 
procedure: 

1. Construct an ordinary allsource, i.e. a unit color-spin vector on every site. 

2. To the allsource, apply the gauge transformation that transforms from the original 
unfixed gauge to the Coulomb gauge. 
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3. Smear the source terms on each site in Coulomb gauge by convoluting with an expo- 
nential smearing function. (This is most efficiently done in momentum space using 
FFT's.) 

4. Transform the smeared allsource back to the original unfixed gauge. 

This smeared allsource can be used as the source for the quark propagator calculation, 
and the subsequent analysis is identical to that of the local hairpin correlator. In Coulomb 
gauge the smeared allsource is a superposition of real exponential sources originating from 
every point. By going back to the unfixed gauge, we attach a random SU(3) gauge phase to 
each exponential, so that a quark loop which starts on one exponential and ends on another 
will have a random phase (even if it actually starts and ends at the same space-time point), 
whereas terms which start and end on the same exponential have no random phase (even 
if they start and end on different points). In this sense, the method is very similar in spirit 



to one introduced earlier by Weingarten, et al |13|, where multiple smeared sources are 



introduced in hadron spectroscopy calculations by attaching random U(l) phases to each 
source. 

The ratio of ground-state overlaps of the smeared and local sources with the pion is 
easily and accurately determined from the behavior of the smeared-smeared and smeared- 
local valence pion propagators. As discussed in the next section, values for the pion mass 
and for the ground-state overlaps are obtained from a combined fit to the propagators using 
smeared pseudoscalar, local pseudoscalar, and local axial vector sources and sinks. 

Define the local and smeared V'T^V' operators by Pi{x) and Ps{x) respectively, and 
measure the corresponding matrix elements, 

{0\P,{0)\7T{p = 0)) = f^^ i = l,s (22) 

To test for the presence of excited states in the hairpin correlator, define the smeared and 
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local hairpin correlators (at zero 3-momentum) A^f^\t), i = s,l and plot the ratio 

If there are no excited states, this ratio should be equal to unity. In FigQ we plot this 
ratio for one of the lightest mass hairpins calculated {k = .1427 or uiq = .0137 in lattice 
units). The absence of any excited state contamination in the hairpin propagator is striking. 
By contrast, the ratio of valence propagators at small times is substantially larger than its 
asymptotic value, indicating that the local valence propagator has a larger excited state 
contribution. We conclude that the hairpin vertex is very nearly decoupled from excited 
pseudoscalar states. In Fig[|lO|, we also plot the results of a similar analysis for a heavier 
quark mass {k = .1410 or ruq = .0559 in lattice units). Here the relative contribution of 
excited states to the valence propagator (also shown in the plot) is even larger than in the 
light mass case. The hairpin propagator, on the other hand, still exhibits little if any excited 
state contribution. For t > 2, there is no significant departure of the hairpin ratio from its 
asymptotic value. This analysis of the smeared-to-local hairpin ratio has been carried out 
at all the other mass values with similar conclusions. In no case is there any significant 
indication of excited states for t > 2. 

The demonstrated absence of excited states from the hairpin diagram allows us to make 
effective use of the time-dependence of the correlator to investigate the structure of the 
hairpin vertex. The simplest assumption, often invoked in phenomenological discussions, 
is that the hairpin is simply a momentum independent mass insertion tjiq. With this 
assumption, the quenched hairpin correlator in momentum space is given by 

^h{p) = fp 2 ^0 2 2 fP (24) 

Fourier transforming over pQ, this implies a time-dependence for the zero momentum prop- 
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Figure 9: Ratio R{t) of the local-source hairpin divided by the smeared- source hairpin 
with K = .1427, normalized by the asymptotic valence propagator ratio. 

agator of 

Ah{p = 0;t) = ^^(1 + m^t)e-^-' + ^ Nt - t) (25) 

This structureless hairpin vertex is suggested by large arguments, but it is important to 
test for the more general possibility that the vertex has some additional dependence. To 
lowest order in a p"^ expansion, this would generalize the above analysis by the replacement 

ml ^ n(p2) ^ U{-ml) + + ml)U'{-ml) + ... (26) 

To test for p'^ dependence of the hairpin vertex, and to estimate its effect on the de- 
termination of the r]' mass, we carried out two sets of correlated fits to the hairpin time- 
dependence, one with the pure double-pole formula ( p5| ) and one to the single-pole -|- double- 
pole formula resulting from ([2^). In all these fits, the pion mass was held fixed at the value 
given by the valence propagator analysis. Since we have already demonstrated that there is 
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Figure 10: Ratio R{t) (boxes) of the local-source hairpin divided by the smeared- 
source hairpin with k — .1410, normalized by the asymptotic valence propagator 
ratio. Also shown ( x 's) is the same ratio for the valence propagator. 
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very little excited state contamination in the hairpin propagator, the range of times used in 
the fits is taken to be 2 < t < 10. To summarize the overall results of these fits, the hairpin 
time dependence for all the quark masses studied is well described by a single-parameter fit 



to the pure double-pole formula (25). (Here the pion mass is not a fit parameter, since it is 
already accurately determined from the valence propagator.) In Fig[|ll] we show an example 
of a pure double pole fit to the hairpin correlator over the entire accessible t range. We 
conclude that the hairpin vertex is reasonably well described by a momentum-independent 
mass insertion. The final results for ttt-q, given in Table 1, are extracted from pure double 
pole fits. Using Equation (||) and the lattice value for /7r(= Ja/V^) (see Table 2), this gives 

5 = .062(7) (27) 

for Csw = 1.57, and 

5 = .044(5) (28) 

for Csw = 0. Note that, in the pure double-pole approximation to the hairpin correlator, 
the value of 5 from xt and that obtained from the hairpin residue are related by a factor 
{fAtn^ /2fpmq)^, which should be unity by the chiral Ward identity. The agreement between 
(p^) and (^) for the clover-improved calculations can be traced to the following two facts: 
(1) The double-pole formula gives a good description of the hairpin correlator for all time 
separations, and (2) The Ward identity is well-satisfied for Csw = 1.57. (By contrast, the 
current algebra quark mass is about 40% smaller than the bare mass for Csw = 0-) 
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Figure 11: One-parameter fit to a pure double- Goldstone pole form for the hairpin 
correlator at k,— .1420. Pion mass is fixed from the valence propagator analysis. 

Table 1: Value of mo, the hairpin contribution to the r]' mass, for Cgw = 1-57 (first 
two columns) and Cgw — (last two columns) at j3 = 5.7, for Nf—3. 





moa 


K 


niQa 


.1410 


.517(23) 






.1415 


.534(24) 


.1630 


.179(10) 


.1420 


.554(24) 


.1650 


.239(16) 


.1423 


.568(26) 


.1667 


.289(13) 


.1425 


.576(27) 


.1675 


.321(13) 


.1427 


.576(30) 


.1680 


.338(13) 


.1428 


.576(33) 


.1685 


.353(18) 




.580(27) 




.393(15) 
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5 Quenched chiral logs in the pseudoscalar mass 

The effect of quenched 77' loops on the chiral behavior of the pseudoscalar mass is one of the 
most definitive predictions of the quenched chiral log analysis. [||, |7| In a previous analysis of 
the pion mass as a function of quark mass,Q] no evidence was found for quenched chiral log 
behavior at /? = 5.7 for unimproved Wilson fermions, with a one-standard-deviation upper 
bound on the chiral log parameter of (5 < .07. This was also shown to be consistent with the 
size of the hairpin propagator. That analysis was done before the development of the MQA 
method for resolving the exceptional configuration problem, and the lightest pion mass used 
was m-j^a = .253 (hopping parameter k = .1680 and Csw = 0). With MQA improvement of 
quark propagators we obtain much better statistical errors on and also are able to go 
to a much lighter quark mass (7717^0 = .164 or k = .1687 for Csw = and rriT^a = .244 or 
K = .1428 for Csw = 1.57). As we discuss in this section, this improved analysis allows us 
to observe clearly the quenched chiral log effect in the pion mass with a value of the chiral 
log parameter 6 slightly less than the previously established upper bound (5 = .054) for the 
Csw = case. The value of 6 is somewhat larger for clover improved quarks {6 = .073), 
suggesting that the suppression of 5 compared to the expected continuum value ~ .17 may 
be at least partially due to finite lattice spacing effects. The recent results from CPPACS of 

6 = .06 — .12[P] for several values of /3 is consistent with this possibility but is not accurate 
enough to observe any clear lattice spacing dependence. 

To extract a value of the chiral log parameter from the pion mass as a function of bare 
quark mass, we shall consider in this section only the case where the two quark masses 
involved are equal (in Section 7 chiral log formulas for unequal quark masses will be derived 
from a model chiral Lagrangian and used to perform global fits to the b lattice data). In 
the equal quark mass case, the chiral logs sum up in a leading log approximation to give an 
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Table 2: Pseudoscalar masses and decay constants for Cgw = 1-57 (first four columns) 
and C«,„ = (last four columns at /? = 5.7. 



n 


mpa 




fpa 


K 


mpa 


fAa 


fpa 


.1400 


.603(2) 


.196(2) 


.458(5) 


.1610 


.647(1) 


.224(3) 


.518(9) 


.1405 


.556(2) 


.190(2) 


.444(5) 


.1630 


.558(1) 


.199(2) 


.471(9) 


.1410 


.505(2) 


.183(2) 


.430(6) 


.1650 


.458(1) 


.174(2) 


.424(8) 


.1415 


.450(3) 


.176(2) 


.418(6) 


.1667 


.356(2) 


.154(2) 


.387(8) 


.1420 


.386(3) 


.169(2) 


.410(7) 


.1675 


.297(2) 


.144(2) 


.371(8) 


.1423 


.342(4) 


.165(3) 


.410(9) 


.1680 


.254(2) 


.136(5) 


.358(12) 


.1425 


.307(4) 


.163(3) 


.413(10) 


.1683 


.221(3) 


.132(6) 


.353(19) 


.1427 


.267(5) 


.161(4) 


.424(14) 


.1685 


.195(4) 


.129(7) 


.342(25) 


.1428 


.245(6) 


.161(5) 


.439(17) 


.1687 


.164(5) 


.126(11) 


.345(38) 






.151(2) 








.122(2) 





anomalous power law dependence of the squared pion mass on the quark mass, 



const. X mq^^ 



(29) 



For Csw = 1-57 we calculated the pion mass at 9 values of hopping parameter ranging 
from K = .1400 to .1428. The masses are obtained from a combined, correlated fit of 
smeared-local and smeared-smeared propagators, using a smeared pseudoscalar source, a 
local pseudoscalar source, and a local axial-vector source. The pion masses obtained are 
listed in Table 2. 

Figures ||12[ and [^] exhibit the chiral log effect in the pseudoscalar mass graphically. 
The first plot includes all nine values of quark mass. The solid line is the best quadratic fit 
(i.e. = Amq + Bm^) to the four heaviest masses (with included as a fit parameter). 



The second plot (Fig|]l3|) is an expanded view of the small mass region. It shows clearly 
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Figure 12: mp vs. quark mass for /3 = 5.7 and Csw = 1-57. Solid line is the second 
order xPT fit to the four heaviest masses. 

that the hght pion masses fah below the quadratic extrapolation of the heavier masses (solid 
line). The dashed line is a fit of the lowest five masses to the chiral log formula (p9|) (again 
with the value of Kc as one of the fit parameters) . We find 

5 = .073(20) (30) 

We have also calculated the pion masses for Csw = on an ensemble of 200 gauge 
configurations at /3 = 5.7 on a 16^ X 32 lattice (the "a" ensemble from the ACPMAPS 
library). Here also we calculate m-,j at nine values of hopping parameter ranging from 
K = .1610 to K = .1687. For this case we were able to go to an even smaller pion mass 
of m-j^a = .164 which is less than 200 MeV in physical units (using the charmonium scale 
a^^ = 1.18 GeV). The pion masses are listed in Table 2. Again the quenched chiral log effect 
is clearly visible, with the lightest-mass points falling significantly below an extrapolated 
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Figure 13: Quenched chiral log fit (dashed hue) to the hghtest 5 mass values with 
S = .073. Solid line is the extrapolated xPT fit to the heaviest four masses. 



quadratic fit. Fitting to tfie quenched cliirai fog formula (29) we find 



6 = .054(20) 



(31) 



Consistent with the direct hairpin caicuiation, the vaiue of the chiral log parameter from 
the analysis is somewhat smaller for Cgw = than for Cgw = 1-57, although in this case 
the error bars are larger, so the difference is only marginally significant. 
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Figure 14: Pseudoscalar decay constant fp as a function of quark mass. Enhance- 
ment at small ruq is a QCL effect. 

6 Chiral behavior of the pseudoscalar and axial- 
vector matrix elements 

As discussed in the Introduction, the chiral behavior of the pseudoscalar and axial-vector 
decay constants fp and /a provide further tests of quenched chiral log predictions. When 
the two quarks in the pseudoscalar meson have equal mass, we should find a clear contrast 
between these two quantities: fp should exhibit a QCL factor oc while should 

have a smooth, nonsingular chiral limit. The values of fp and are obtained from the 
combined fit to smeared-smeared and smeared-local propagators discussed in Section 5. 
The effects of full tadpole improvement for the normalization of quark masses and quark 
operators have been included in our calculations, but we have not carried out a complete 
nonperturbative 0(a) improvement program. The numerical results are presented in Table 
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Figure 15: Axial-vector decay constant as a function of quark mass. Note the 
absence of a QCL enhancement, as predicted by theoretical arguments. 

2. For the Csw = 1-57 results, the mass dependence of the decay constants is shown in 
Figures and ||l^. It is clear from these plots that at the very lightest masses, the value 
of fp is significantly larger than a linear extrapolation of the heavier mass results, consistent 
with the singular (m^)^'^ expected from QCL effects, while /a exhibits no sign of singular 
behavior, and is well described by a linear fit. In both figures, the solid line is the best 
linear fit to the four heaviest masses. 

Since the singular behavior of fp in the chiral limit is determined by the ubiquitous 
chiral log parameter 5, it should be possible to use the results for fp to obtain another 
estimate of this parameter. The extraction of a value of 6 is made somewhat more difficult 
by the fact that, in addition to the singular QCL effect, the mass-dependence of fp also 
has a significant contribution from higher order chiral perturbation theory (i.e. terms linear 
in m^). With the accuracy of our data, the two contributions are difficult to disentangle. 
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Figure 16: Ratio fpj Ja as a function of quark mass. Solid line is a QCL fit with 
b = .085. 

Although the QCL effect is clearly visible, a fit to fp which includes both QCL terms and 
xPT terms is rather unstable and the resulting value for 6 is poorly determined. We can 
do much better if we make an additional, phenomenologically motivated assumption that 
the perturbative slopes of fp and /a are approximately equal 10. The data in Figures 
1 14] and are consistent with this assumption. The ratio fp/ Ja should thus exhibit a 
relatively pure chiral log behavior. 



fp 

— = const. X {m^)~^ 
JA 



(32) 



The lattice results for this ratio are shown in Fig|16|, along with the best fit to the QCL 
formula (^). This gives a value of the chiral log parameter of 



5 = .085(23) 



(33) 
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7 Extraction of 6 from mass and decay constant 
cross-ratios 



To facilitate a comparison with previous work by the CPPACS cohaboration we have 
extracted the chiral log parameter 5 from our full set of b-lattice clover-improved results 
for masses and decay constants of the 45 independent mesons which can be formed from 
the nine available quark masses, using the cross-ratio method introduced in Q. For a given 
meson parameter yij (here i,j label the quarks in the meson and run from 1 to 9), the 



cross-ratio Rij is defined as follows 



Rij ^ (34) 



Let Mij,fp-ij, and fA;ij denote the mass, pseudoscalar and axial- vector decay constant 
of the meson with quark content i,j (and quark masses rrii and nij). Then, with either 
Hij = Mfj/{mi + mj) or y^j = jj^, one has, to leading order in J, hut ignoring higher order 
chiral perturbation theory contributions (this amounts to setting L^,Ls = in the chiral 
formulas ( ^^J59| ) discussed in Section 8 below), 

Rij = const. X {1 + 6Xij) (35) 

where the chiral logarithm is contained in the factor Xij. At infinite volume, this factor 
becomes 

Xi,- = 2 + ^ In (36) 

TUi — ruj rrii 

Our lattices are at smaller volume that those in the work cited above (physical extent 
~ 2 fermi as compared to 3 fermi in j^) and we go to considerably smaller quark masses, 
so we have used a finite volume version of the fitting parameter Xij (see Section 8 below 
for a discussion): 

Xij = '^lij In Ijj (37) 
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Figure 17: Cross- Ratios of Mf^/ {mi + mj) as a function of finite volume Xij variable 

where the finite volume sums lij are defined in (^^. 

Fitting the cross-ratios of the masses of all 36 off-diagonal mesons, we obtain 

<5 = 0.0604 ± 0.0155 (x^ = 54/35 d.o.f ) (38) 

while a fit of cross-ratios of the decay constant ratios gives 

J A;ij 

6 = 0.0708 ± 0.0104 (x^ = 37/35 d.o.f) (39) 

Finally, a combined fit of both the mass and decay-constant ratios gives 

6 = 0.0732 ± 0.0081 (x^ = 118/71 d.o.f) (40) 

The expected linearity in the X-variable of ( |35| ) is displayed graphically in Fig|l^ for the 
pseudoscalar masses and in Fig[^ for the ratio of pseudoscalar to axial decay constants. 
The results of fully correlated fits are displayed as solid lines. 
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Figure 18: Cross- Ratios of fp-ij/ fA;ij as a function of finite volume Xij variable 
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8 Comparison with quenched chiral perturbation 
theory 



For pions made from a quark and antiquark with unequal masses, the form of the quenched 
chiral log effect is more complicated ||3|, |^. For the range of masses we consider, it is sufficient 
to keep only lowest order terms in 6 (i.e. one- loop terms), or equivalently, in a hairpin mass 
term which can be included explicitly as a correction to the basic chiral Lagrangian : 

2 



where 



+ LsTT{x^Ux^U + U^xU^x)+Ci,^irpin (41) 



/:hairpin = -^ml^{iTvln{U^) - iTTln{U)f (42) 
z o 



The lowest order chiral Lagrangian has been supplemented by the chiral symmetry breaking 
terms -L5 of O(p'^m^) and Lg of O(m^) which model the leading mass-dependence of the 
slope in the pseudoscalar masses. Starting from this Lagrangian, we can derive explicit 
formulas for the pseudoscalar masses, and pseudoscalar and axial vector decay constants, 
consistent through order and including the effects of the hairpin mass insertion (assumed 
local) through the term £hairpin in (|4l]). The coefficients L^,Ls follow the notation of 



Gasser and Leutwyler |14]. The evaluation of the one loop chiral integrals appearing in 
this calculation has also been carried out at finite volume (appropriate for the physical size 
of the lattices used), and with a Pauli-Villars subtraction to regulate the ultraviolet. For 
example, logarithmically divergent integrals such as 



If 1 1 

^iJ = ^J d^Pp2 + Mf p2 + ^^^^ 
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are replaced by 

= 167r2 Y^iDip, Mi)Dip, Mj) - D{p, Af) (44) 

p 

where the momentum integration is now a discrete sum over a finite volume free boson 
propagator D{p, M). We have typically chosen the cutoff scale A ~ i but sensitivity of the 
results to this choice is very small. One also encounters quadratically divergent graphs in 
the course of the calculation, which are regulated as follows: 

h = 167r2 Y.{D{p, Mi) - D{p, A) - (A^ - Mi)D{p, Af) (45) 
p 

With these preliminaries, we find the following expression for the pseudoscalar masses 
(squared), up to first order in the hairpin mass and (independently) in L5 and Lg: 



Mij = ^4^(1 + W.,)x 



1 

1 

'P 

Jij = {Ii + Ij-{Mi + Mfj)Iij)/2 (47) 



{1 + -^(SLg - 4L5)(xi + Xj) + pSiSLs - 4.L^){2IuXi + '^hiXj + {Xi + Xj)Iij) 
+ S^L^SJij} (46) 



The quantities Xi encode the quark masses: our data includes values for 9 different kappa 
values, so the indices i,j above run from 1 to 9, allowing for 45 independent quark-antiquark 
combinations. Thus 

Xi = 2romi (48) 

where tq is a slope parameter (which we also extract from the fits) and we have used the 
pole value for the quark mass: 



^ ^^^^ + 27 " ^^^^ 
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Similar formulas were obtained for the pseudoscalar and axial decay constants fp-ij , fA;ij 
and are listed in Appendix 2. We have performed global fits to the masses and decay con- 
stants for all 45 mesons in order to extract the parameters f,rQ,6,L^ and Lg. The fits 
were performed for a variety of time-windows (for the 12^x24 b lattices, on time windows 
5-11,6-11,7-11 and 8-11) in order to isolate any remaining sensitivity to higher state con- 
tamination, and the results for the various chiral parameters, as a function of the initial 



time for the window, are shown in Fig|19]. With only 300 independent configurations, it is 
not possible to obtain a sufficiently stable covariance matrix to fit all 135 masses and decay 
constants, so these results reflect an uncorrelated fit to all meson parameters using (p6|), 
(H) and (H). 

To summarize our results, a global fit to all pseudoscalar masses and decay constants, 
using a time window of 6-11, gives a final value for the chiral log parameter of 

5 = 0.059 ± 0.015 (50) 

while the slope and critical kappa parameters in (48,49) are determined as 



ro = 1.99 ±0.12 (51) 

Kc = 0.143287 ± 0.000022 (52) 
For the chiral breaking parameters L5 and Lg, our fits give 

L5 = (2.5 ± 0.5) X 10"^ (53) 

and 

Lg- ^Ls = (0.2 ±0.1) X 10-3 (54) 

The dimensionless chiral parameters L5, Lg are only roughly determined by phenomenology. 
Recent estimates [^] give 

L5{Mp) = (1.4 ± 0.5) X 10-3 (55) 
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Figure 19: Global Fits, using (^ , |58| , |59D , for Kc, S,rQ, L^^g, /a = \/2/,and fp = r^f / 
Results are shown for time windows 5-11,6-11,7-11, and 8-11 
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renormalized at the rho mass, with the combination Lg — consistent with zero: 

1 

Lg - -L5 = (0.2 ± 0.4) X 10-2 (56) 

Finally, our result for the axial decay constant corresponds to a value for the pion decay 
constant (bare, in lattice units) of 

f = U = 0.1066 ± 0.0024 (57) 
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9 Summary and discussion 

The calculations presented above confirm all the essential features of anomalous quenched 
chiral behavior suggested by continuum calculations. Our lattice studies lead to the follow- 
ing basic conclusions: 

1. The MQA pole-shifting technique allows for accurate computation of meson and hair- 
pin correlators down to small quark masses (iW^ < 200 MeV). Probing this region is 
essential in order to obtain reliable signatures of anomalous chiral behavior. 

2. The hairpin vertex has only a small coupling to excited states, and very gentle mo- 
mentum dependence. This suggests that it may be accurately modelled by a local 
mass-insertion term in a chiral Lagrangian. 

3. Determination of the chiral log parameter 5 using five separate methods gives con- 
sistent results. A summary of our results for this parameter indicating the various 
methods used is displayed in Table 3. The overall average of our Cgw = 1-57 methods 
gives S = 0.065 it 0.013. Our values for 6 at (3 = 5.7 are considerably smaller than 
those expected from a naive continuum analysis, but are in agreement with previous 
lattice estimates. 

4. Using the MQA technique, meson properties (masses and decay constants) can be 
extracted with sufficient accuracy to allow a fit of higher order chiral parameters, 
such as L5 and Lg. 

Careful quantitative studies of chiral behavior in quenched QCD in comparison to 
quenched chiral perturbation theory can provide a great deal of insight into the connection 
between QCD and the effective chiral Lagrangian that describes its long range behavior in 
the limit of small quark mass. Even if the numerical simulation of full QCD were not so 
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Table 3: Summary of results for chiral log parameter 5 



Method 


C 


5 


Witten-Veneziano 


1.57 


0.063(6) 


liMiri")in \^('r(.(^x 


1.57 


0.062(7) 


diagonal mesons 


1.57 


0.073(20) 


ratio fit, masses 


1.57 


0.060(16) 


ratio fit, fp/fA 


1.57 


0.071(10) 


global fit 


1.57 


0.059(15) 


Witten-Veneziano 





0.074(9) 


hairpin vertex 





0.044(5) 


diagonal mesons 





0.054(20) 



expensive computationally, the study of chiral behavior in quenched QCD would still be 
of theoretical interest. For example, the Witten-Veneziano relation connects the r/' mass 
to the topological susceptibility of quenched QCD. The geometric summation of multiple 
•q' mass insertions is only the simplest example of how, in some cases, the most important 
effects of the full QCD fermion determinant can be incorporated into a quenched result, 
with the guidance of chiral perturbation theory. Because of the smallness of the parameter 
5, QCL effects are adequately described by one-loop xPT, even for quark masses close to the 
physical up and down mass. It is thus straightforward to apply appropriate and calculable 
QCL corrections to quenched results. Of course the masses, decay constants, and higher 
order chiral Lagrangian coefficients obtained in quenched QCD will differ somewhat from 
those of the full theory, but all of the disturbing structural properties of the quenched the- 
ory (lack of unitarity, absence of topological screening, etc.) can be systematically repaired 
in the context of chiral perturbation theory. Further precision studies of anomalous chiral 
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behavior in the quenched meson and baryon spectrum should provide additional insight into 
the origin and structure of chiral symmetry in QCD. The results presented in this paper 
provide strong support for the usefulness of the MQA technique to facilitate these studies. 

Ideally, similar studies should be performed using an exactly chirally symmetric Dirac 
operator which satisfies Ginsparg- Wilson relations, e.g. the Neuberger operator |jl^. Ex- 
plicit Ginsparg- Wilson chiral symmetry would resolve the exceptional configuration problem 



ab initio. Unfortunately, such operators are necessarily not ultralocal|16], and are difficult 
to invert or diagonalize numerically. 

The MQA method attempts to account for the most salutary effect of an explicitly 
chirally symmetric approach. The underlying assumption in this procedure is that the 
most important effect of the explicit chiral symmetry breaking contained in the Wilson- 
Dirac operator is the real displacement of its small eigenvalues. The MQA method merely 
removes these displacements in a compensated manner. The principal disadvantage of the 
method is its apparent lack of locality. Using a basis of hopping terms (nearest-neighbor, 
next-nearest neighbor, etc.), it may be possible to identify terms in an ultralocal expansion 
of the Dirac operator which correspond most closely to an MQA improved Wilson-Dirac 
operator. The additional hopping terms would have the effect of reducing the dispersion 
of the small real eigenvalues as well as inducing small rotations of the basis wavefunctions, 
etc. Such an analysis has not been carried out yet. The clear quantitative success of the 
MQA procedure in restoring desired chiral behavior is a promising indication that such 
eigenmode-based methods can be both efficient and effective in removing the dominant 
spurious effects of chiral symmetry breaking contained in the Wilson-Dirac formulation of 
lattice fermions. 
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Appendix 1: The allsource method 



The allsource method, first applied to the rj' mass calculation in Ref. |11], is used here 
to calculate hairpin diagrams as well as the pseudoscalar charge for the determination of 
winding numbers. We have also introduced a generalization of this allsource method which 
allows the calculation of closed loops which originate from a smeared source, i.e. the two 
ends of the quark propagator are contracted over color and spin, but are spatially separated 
with an exponential weight function. The method relies on the fact that gauge noninvariant 
terms in the calculation will cancel out due to random SU(3) gauge phases. For example, 
in the calculation of a single quark loop, the closed loop terms where the quark starts and 
ends on the same point (or on the same exponential source in the smeared calculation) add 
coherently when summed over sites, while loops which start and end on different sources 
have random phases and cancel. 

While the cancellation of random gauge phases in the allsource method should work 
arbitrarily well for a large enough ensemble of gauge configurations, it is very instructive 
to test this method in a situation where we know the exact gauge invariant answer which 
can be used to check the accuracy of the random phase cancellation. To carry out such a 
comparison, we have calculated the ordinary valence pion propagator using the allsource 
technique. This calculation uses the same allsource propagators that were used to calculate 
hairpin correlators. In the hairpin calculation, we computed the correlator of two closed 
loops by contracting the two color indices of each allsource propagator with a fixed time 
separation between the two propagators. To calculate the valence propagator, we instead 
take the same two allsource propagators with fixed time separation and cross-contract the 
color indices to form a single loop from the two propagators and then project out the color 
singlet component. (The last step amounts to leaving out the terms in which all four color 



indices are equal and then multiplying by a factor of 3/2.) In Fig[20|, the valence pion 
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Figure 20: Comparison of valence pion propagator computed by allsource method 
(x's) and conventional method (boxes). Data points are offset slightly for clarity. 

propagator calculated by the allsource method is compared with the results of the standard 
calculation using local and smeared sources on a fixed timeslice. The results are shown for 
hopping parameter k = .1425, but similar agreement is obtained at all kappa values. 

The agreement is excellent and well within statistical errors. For small times {t < 5), 
the errors on the allsource calculation are actually smaller than those from a fixed source. 
(Remember that the allsource calculation allows the meson propagator to be averaged over 
all locations on the lattice, thus increasing the effective statistics relative to the fixed source 
method.) Unlike the fixed source calculation however, the statistical errors for the allsource 
calcTilation arc morc-or-less constant in absolute magnitude for all time separations. This 
results in a signal-to-noise ratio that gets rapidly worse as we go out in time, just as in the 
hairpin calculation. This roughly constant noise level is presumably the effect of incomplete 
cancellation of random gauge phases. Because of this, the allsource method is an inferior way 
of studying the asymptotic behavior of the pion propagator. Nevertheless, for short times 
it accurately reproduces the results of the standard method. It is a fortuitous circumstance 
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that the hairpin correlator is found to be almost entirely free of excited state contamination. 
This allows us to extract the ground state vertex insertion from its value at relatively short 
times where the allsource technique is accurate. 
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Appendix 2: Quenched chiral results for decay 
constants 

A next to leading order chiral perturbation theory calculation of the pseudoscalar and axial 
decay constants can be carried out along the same lines as discussed in Section 8 for the 
pseudoscalar mass spectrum. Starting from the Lagrangian and using the notations 
introduced there, we find 

/P;ij = \/2/ro(l + j6{Iii + Ijj + 2Iij)) X 



{ 1 + -^(SLg - 2L5)te + Xj)(l + <5(2Iia* + + te + Xj)Iij)/{Xi + Xj)) 

2 

-j^L^6{IiiXi + IjjXj + iXi + Xj)Iij) 
+ j2^5S{Jii + Jn)} (58) 



for the pseudoscalar decay constant, and 

fA;ij = V2f{l + ^6{Iii+I,j-2Iij))> 



{ l + j^L5{xi + Xj) + 2j^L56{2{xi + Xj)Iij-XiIii-XjIjj)} (59) 



for the axial decay constant. These formulas, together with 



M^. = ^^(1 + W.,)x 



{1 + (BLg - 4L5){x^ + Xj) + j^^iSLs - AL^){2IuXt + '^IjjXj + {x^ + Xi)^*j) 
+ S^L^SJi,} (60) 



for the pseudoscalar masses, were used to perform the global fits described in Section 
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